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ABSTRACT 



This thesis describes an evaluation of the active and 
passive acoustic detection modules of the Naval War College 
Gaming System (NWCGS) that has been installed at the Naval 
War College in Newport, Rhode Island. The specific intent 
of the evaluation is to verify that the model is theoretically 
sound. This evaluation compares the NWCGS model to other 
existing acoustic detection models. Recommendations for 
improvement to the model are also presented. 
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THE NAVAL WAR COLLEGE GAMING SYSTEM MODEL 



The detection modules used by the Naval War College 
Gaming System (NWCGS) simulate the acoustic detection proc- 
ess using the generalized expected value sonar equations 
developed during the Second World War. Signal Excess, the 
difference between the existing signal-to-noise ratio and 
that required to achieve a specified level of performance, 
is the fundamental parameter that determines whether or not 
a detection has occurred. In the NWCGS , the value of signal 
excess is determined by the following two equations: 

SE = SL - TL + DI - NL - DT (li 

when modeling mpassive sonar detections, and 

SE = SL - 2TL + TS - (NL + DT) (2) 

when modeling active sonar detections. 

The parameters referenced by equations (1) and (2) are- 
defined as follows: 

SE - Signal Excess 

SL - The source level of the target as measured one 
yard from the acoustic center of the effective 
sound source 

DI - The directivity index of the hydrophone 

TL - The one-way transmission loss from the target to 
detector 
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NL - The non-target noise in the immediate area of 
the receiver 

TS - The target strength as measured one yard from 
the acoustic center of the target 

DT - The detection threshold. This is the signal- 
to-noise ratio (in decibels) required to 
achieve a specified probability of detection 
at a specified probability of false alarm 

A more complete description of these parameters is presented 

in Appendix A. The discussion which follows describes the 

methods used in the Naval War College Gaming Systran to 

determine appropriate values for the parameters of. equations 

(1) and (2) . 

A . METHOD 

The primary method of data retrieval used r>y the NNCG 
system is the multi-parameter "look-up" table. A look-up 
table can be described as an array having indices that 
reference its entries. The entries in the look-up table 
describing acoustic transmission loss are transmission loss 
values for different ranges and different environmental 
conditions. It is an example of a two-dimens icnal look-up 
table. The first index of the table identifies tie environ- 
mental conditions assigned to the area of the game. The 
second index identifies the distance traveled by sound waves 
from the target to the detector. The advantage of a look-up 
table is its ability to make use of calculations done in 
advance. By using a look-up table, the amount of computations 
done during the course of a game is reduced, and the maximum 



7 



possible game speed is increased. The benefits of a faster 
game are not without some cost. The expense of using a look- 
up table is the requirement for additional memory and 
associated hardware used to store the tables. Because a 
look-up table is a discrete approximation of a continuous 
process, the accuracy of any signal excess model that uti- 
lizes a look-up table is limited by the resolution of the 
tables being used. In the transmission less look-up table 
example, the look-up table index describing tc.rget-to-sensor 
range has a resolution of 1000 yards. This means that 
ranges are rounded to the nearest 1000 yards before entering 
the transmission loss look-up table. Increasing the table 
range resolution from 1000 to 500 yards requires twice the 
memory currently required. The presentation that follows 
examines each parameter of equations (1) and [ 2 ), and the 
tables used by the NWCGS to determine values for each. 

B . TRANSMISSION LOSS 

Transmission loss (described in tne NWCGS documentation 
as propagation loss PL) is determined through a two-parameter 
look-up table. The first parameter determines the environ- 
mental conditions in the area of the game. All players are 
assumed to be contained in a single environmental area, and 
the sound velocity profile is assumed to be the same through- 
out all parts of that area. The second parameter of the table 
is range from target to source (rounded to the nearest 1000 
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yards) . For a given environmental area, transmission loss 
becomes a function only of range. Such a table can adequately 
account for acoustic spreading losses, however as described 
in Appendix A, absorption and ocher forms of transmission 
loss depend upon other factors in addition to range. 

1. Direct Path Mode 

When both the target and the searcher are in the 
mixed surface layer, transmission loss is primarily a func- 
tion of range r from target to searcher, target depth D fc , 
and layer depth D 1 . Sound spreads spherically until it fills 
the layer (at transition range r ) , and then spreads cylin- 
drically. Acoustic waves can be scattered at the surface, 
or may leak out of the bottom of the nixed layer. The total 
loss due to surface reflections is equal to the loss per 
surface reflection multiplied by che number of surface 
reflections that occur between car get and searcher. The 
number of reflections occurring is sinply the distance between 
target and detector divided by the distance between reflec- 
tions. The complete expression for direct path transmission 
loss (when both target and searcher are within the mixed 
surface layer) is given by the expression [Ref. 1]: 

TL = 10 log r t + 10 log r + ar + br/r r (3) 

where : 



r = horizontal range from target to searcher 
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r fc = 115 sqrt ID 2 / (D t ~D 1 ) 

a. = absorption coefficient 

b = loss (dB) /reflection = 0 . 63*freq* (1 . 4 ) (sea state) 
= 919 sqrt^D^) 

The transmission loss look-up table used by the NWCGS model 
describes only losses due to cylindrical spreading. Values 
obtained from such a table can only be used to determine values 
for the first two terms of equation (3) . The table does not 
depend upon the environmental area of the ocean (since cylin- 
drical spreading is the same in all environments) ; and it 
can be reduced to a one-dime isional table dependent only 
upon the target-to-searcher range. The curve of Figure 1 shows 
how cylindrical spreading loss varies as a function of range. 

This plot shows than a look-up table having a 1000 yard reso- 
lution can be responsible for table round-off errors that 
exceed 4 dB. In some cases, an error of this magnitude may 
not be acceptable. 

The two transmission loss terms of equation (3) that 
cannot be described by the two-dimensional look-up table 
used in the NWCGS are the losses due to absorption and the 
losses which occur as a result of scattering of acoustic 
energy at the sea surface or "leakage" of energy out of the 
bottom of the layer. Scattering losses are functions of 
range, sea state and signal frequency. A three-dimensional 
look-up table that would be necessary to describe scattering 
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Figure 1. Spreading Loss Profile 



losses is not an efficient use of available storage, so 
execution time calculation of this loss factor is recommended. 
Some execution time speed can be realized by pre-computing 
the parameter 

b = 0 . 63*freq (1.4) (sea state) (4) 

The sea state of a given area remains constant in a given 
area, so b can be precomputed for a number of frequencies 
and stored in a one-dimensional look-up table that is asso- 
ciated with an environmental area. Table 1 is an example 
of such a table using a sea state of 4 . 

2 . Absorption Loss 

To account for the additional transmission loss due 
to absorption of sound by the medium, an absorption loss 
factor should be included. The absorption loss corrections 
currently used in the NWCGS model is applied only to the 
active model. The correction used is of the form: 

TL = TL + 0.033 * f req'*' * * range (Kyd) (5) 

This equation, cited in Reference 2 is derived from research 
done by M.J. Sheehy and R. Halley [Ref. 3] in 1956. Their 
research attempted to determine the attenuation coefficient 
for low frequency sounds . The experiments involved detonation 
of a nuclear device several hundred miles off the coast of 
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TABLE 1 



(sea state 
b 

2.42021 

4.84042 

7.26062 

9.68083 

12.101 

14.5212 

16.9415 

19.3617 

21.7819 

24.2021 

26.6223 

29.0425 

31.4627 

33.8829 

36.3031 

38.7233 

41.1435 

43.5637 

45.984 

48.4042 

50.8244 

53.2446 

55.6648 

58.085 

60.5052 

62.9254 

65.3456 

67.7658 

70.186 

72.6062 



= 4) 




freq 


b 


31 


75.0264 


32 


77.4467 


33 


79 . 8669 


34 


82.2871 


35 


84 . 7073 


36 


87.1275 


37 


89.5477 


38 


91.9679 


39 


94 . 3881 


40 


96 . 8083 


41 


99.2285 


42 


101.649 


43 


104.069 


44 


106.489 


45 


108.909 


46 


111.33 


47 


113.75 


48 


116.17 

0 


49 


118.59 


50 


121.01 


51 


123.431 


52 


125.851 


53 


128.271 


54 


130.691 


55 


133.111 


56 


135.532 


57 


137.952 


58 


140.372 


59 


142.792 


60 


145.212 
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southern California, and recording the acoustic intensity of 
frequencies between 20 and 200 hertz. The data, collected 
from several locations in the eastern Pacific Ocean, was 
fitted using a least-squares regression technique to an 
approximating function described by equation (5). Extrapo- 
lating equation (5) to frequencies above 200 hertz and com- 
paring the predicted results to measured data, their report 
concluded that equation (5) was a reasonable approximation 
for frequencies up to 60 Kilohertz. 

Subsequent studies, cited by Urick [Ref. 4] have 
developed somewhat different approximating functions for 
the absorption coefficient. The most current of these is 
given by: 

-2 ^T^ -2 f^ 

a = 1.86 x 10 T + 2.86 x 10 dB/Kyd (6) 

f 2 + f 2 f t 

Figure 2 is a plot of absorption loss as a function of range 
The two curves depicted are drawn using two approximating 
functions (equations 5 and 6) . Figures 3 and 4 show how the 
choice of approximating function can affect the total absorp 
tion loss correction factor for a given range and frequency. 
It is evident from these plots that the differences between 
the two approximating functions are not large, however as 
the range or frequency increases, the difference between the 
two increases. Assuming that the approximating function 
cited by Urick is based on a larger set of experimental data 
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Figure 2. Absorption Loss Profiles 
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Figure 3. KWCGS Absorption Loss Characteristics, Curve 
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Figure 4. Urick's Absorption. Loss Characteristics Curve 



it is recommended that equation (6) replace equation (5) 
in the NWCGS detection modules. For the sake of consistency, 
it is further recommended that the absorption loss correc- 
tion be applied to acoustic signals in the passive acoustic 
detection modules as well as the active acoustic detection 
modules. In the case of active detections, the frequency 
"f" of equation (6) is the ping frequency of the sonar trans- 
ducer. In the case of passive detections, ”f" should be the 
center frequency of the bandwidth in which the searcher is 
looking. Since equation (6) is a function of frequency 
alone, a one-dimensional absorption coefficient look-up table 
can be constructed. Table 2, derived from equation (6), 
provides data that may be useful in the construction of such 
a look-up table . 

3 . Convergence Zone Mode 

If convergence zone transmission paths are possible, 
the NWCG system is able to adjust the transmission loss to 
reflect this type of contact. Whenever the searcher is in 
the target's convergence zone area, the NWCG system multi- 
plies the unadjusted transmission loss by a convergence zone 
correction factor. In order to reduce the unadjusted trans- 
mission loss to reflect a convergence zone gain, the multi- 
plicative factor must be a real number between zero and one. 
This method is significantly different than the method used 
by Coppens [Ref. 5] to describe this gain (see Appendix A). 
Since this method of adding rather than multiplying the 
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TABLE 2 



(salinity = 35) 
(temperature = 50 C) 



alpha 


freq 


alpha 


.644164 


31 


15.0205 


1.27433 


32 


15.378 


1.89095 


33 


15.7298 


2.49446 


34 


16.0758 


3.08527 


35 


16.4163 


3.66378 


36 


16.7514 


4.23037 


37 


17.0813 


4.7854 


38 


17.406 


5.32922 


39 


17.7257 


5.86218 


40 


18.0404 


6.38459 


4.1 


18.3504 


6.89676 


42 


18.6556 


7.39899 


4 3 


18.9563 


7.89157 


44 


19 .2525 


8.37477 


45 


19.5443 


8.84886 


46 


19.8318 


9.31409 


47 


20.1152 


9.77072 


48 


20.3944 


10.219 


49 


20.6696 


10.6591 


50 


20.9409 


11.0913 


51 


21.2083 


11.5157 


52 


21.472 


11.9327 


53 


21.732 


12.3424 


54 


21.9884 


12.7449 


55 


22.2413 


13.1405 


56 


22.4907 


13.5294 


57 


22.7367 


13.9116 


58 


22.9794 


14.2875 


59 


23.2188 


14.657 


60 


23.4551 
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correction factor has become an accepted standard, it is 
recommended that the NWCG system be modified to take advan- 
tage of the availability of existing data bases. 

4 . Bottom Bounce Mode 

Transmission loss is adjusted in the NWCG system 
whenever the environment supports bottom bounce. A bottom 
bounce flag is set if the range from target to searcher is 
within a bottom bounce region, and all transmission loss 
calculations are multiplied by a bottom bounce correction 
factor. The areas where bottom bounce detections are con- 
sidered likely are those between the maximum direct path 
range and the minimum convergence zone range. While bottom 
bounce transmission paths can extend into both direct path 
and convergence zone areas, the transmission loss for bottom 
bounce is generally larger than tha.t of the other two paths. 
Consequently, any detection that occurs in one of the other 
transmission areas will most likely not be due to bottom 
bounce. The current literature discussed in Appendix A 
describes the following procedure to adjust the transmission 
loss : 

First, increase the effective range from target to 
searcher by the equation: 

R = R/ ( cos 0) (7) 

where 0 is the angle made between the incident sound wave 
and the bottom. The increase in range is due to the indirect 
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transmission path. Second, increase the total transmission 
loss by adding a constant factor BL due to absorption and 
scattering of the signal by the ocean bottom. 

The spreading loss associated with bottom bounce 
transmission paths is spherical. Since the NWCGS spreading 
loss table describes cylindrical spreading, spreading loss 
values obtained from this table must be doubled to reflect 
that difference. The loss associated with a bottom reflec- 
tion is an environmental parameter that should be included 
in the table that describes the environment. 

An important aspect of bottom bounce often overlooked 
is its value in cases of layer depth interdiction. When the 
searcher is above the mixed layer, and the target is below 
that layer (or vice versa) direct path ranges are severely 
limited due to the refracting effects of the layer. Bottom 
bounce signals have much higher angles of incidence, and are 
therefore less affected by refraction. If layer depth inter- 
diction occurs in bottom bounce areas, detections should be 
based strictly on bottom bounce transmission paths. 

5 . Conclusion 

Figure 5 shows typical transmission loss profile for 
the three transmission modes (direct path, bottom bounce, 
and convergence zone) . Using the equations of the NWCG sys- 
tem, a similar transmission loss profile has been computed 
and shown in Figure 6. The two figures are significantly 
different, especially in terms of bottom bounce transmission 
loss. A more extensive transmission loss look-up table is 
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Figure 5. Transmission Loss Profile 
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Figure 6. NWCGS Transmission Loss Profile 



required to adequately describe spreading, absorption, trans- 
mission path, etc. This table could eliminate the need for 
bottom bounce and convergence zone correction factors , and 
could be used to describe environmental phenomena that are 
difficult to describe analytically. The NWCGS transmission 
Loss look-up table could be modified to include frequency, 
in addition to environmental area and range, as an entering 
argument to the table. Additionally, the resolution of the 
table (currently 1000 yards) could be increased. 

If the current structure of the transmission loss 
look-up table cannot be altered, the following recommenda- 
tions are made: 

1. Use the propagation loss table currently in the NWCG 
system to describe transmission loss due to cylindri- 
cal sreading only. Spreading losses are a function 
of range only, and the difference between cylindrical 
and spherical spreading losses is a multiplicative 
factor of two. All transmission paths are subject 

to spreading losses, and therefore all transmission 
paths will reference this table. 

2. For direct path transmissions, transmission loss should 
be described by the equation: 

TL = 10 log + 10 log r + ar + br/r r (3) 



where: 



r 



t 



115 sqrt (D^/(D t -D tl ) 



r 



r 



918 sqrt(D t ) 



a 



absorption loss/yard 



b 



loss (dB) /ref lection 



0.63*(1.4) 



(sea state) 



* (freq) 



The value r^ is fixed by the environment selected. 
Consequently, it is possible to compute the value of 
10 log (r t ) beforehand, and include it in the parameter 
list of the environmental area. The value of 10 log (r) 
is obtained directly from the propagation loss look-up 
table. Absorption loss per yard "a" can be determined 
by a call to another look-up table (see Appendix A) . 

The loss per bounce parameter "b" is a function of both 
sea state and frequency; and a two-dimensional look-up 
table similar to Table 2 can be constructed and used 
to determine "b". 

3. For bottom bounce transmission paths, the transmission 
loss should be reflected by the equation: 



The range from target to searcher is divided by cos 9 
to account for the increased transmission path. The 
adjusted range is an entering argument for the propagation 



TL = 20 log (r/cos 0) + ar/cos 0 + BB 



(9) 
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loss look-up table. Since bottom bounce spreading 
is spherical rather than cylindrical, the value 
obtained from the table should be doubled. Absorption 
coefficient "a" is the same as previously described, 
and the bottom bounce loss BB is a parameter determined 
by the environmental area. 

4. In areas wihtout bottom bounce, if the detector is 

above the mixed layer and the target is below (or vice 
versa) , it is reasonable to assume no detections are 
possible. Under the same circumstances in a bottom 
bounce area, detections via bottom bounce transmission 
paths are possible, and should be investigated. This 
would require the modification of the layer depth 
interdiction test. If the detector is above the layer 
and the target is below (or vice versa) and the bottom 
bounce available flag is not set, then no detection 
can occur. If the bottom bounce available flag is 
set, transmission loss is calculated for bottom bounce 
mode . 

C. THRESHOLD CROSSING MODELS 

A major use of the passive or active sonar equation is 
to determine mean signal excess. Mean signal excess is the 
average difference (in decibels) between the measured signal- 
to-noise ratio and the signal-to-noise ratio which will yield 
a 50 percent probability of detection for a specified proba- 
bility of false alarm. The NWCGS utilizes a signal excess 
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threshold crossing model to account for the randomness in 
the components of the sonar equation, and to determine the 
probability of detection. This model defines the signal 
excess at a time t to be a random variable X SE (t) and a 
detection to be the event {x gE (t) >_ 0}. The signal excess 
at time t can be expressed as : 



where SE(t) is the mean signal excess, and X(t) is a random 
variable that determines the variability of the signal 
excess at time t. Values for X(t) are determined by drawing 
from a normal distribution having a mean of zero, and a 
variance that is determined by the variability of the com- 
ponents of the sonar equation. The probability that a 
detection will occur can be written as follows: 



Since X ( t) is normally distributed, this implies that 



The mean signal excess is used as an entering argument for 
an inverse normal look-up table that determines the proba- 
bility of detection. 




SE(t) + X(t) 



( 11 ) 



P(X SE (t) >0} = P{ X ( t) > -SE(t)} 



P{X SE (t) >0} = ${SE ( t) /a} 



(13) 
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Two modes are available for converting the detection 
probability into a detection decision, the deterministic 
mode and the probabilistic mode. In the deterministic mode:, 
the probability of detection is compared to a probability 
of detection threshold that has been set by the umpire during 
the initialization of the game. If the computed probability 
of detection exceeds this predetermined threshold, a detection 
is recorded. The alternative to the deterministic mode is 
the probabilistic mode. In this case, the probability of 
detection is computed, and the detection process is simu- 
lated by a draw from a Uniform (0,1) distribution. A detec- 
tion occurs whenever the random variable drawn is less than 
the computed probability of detection. 

When using the deterministic mode with a given standard 
deviation a, the probability of detection threshold initialized 
by the game director defines a signal excess to sigma ratio 
threshold (SE/a) TH by: 

$(SE/a) TH = (Pd) TH (14) 

If the probability of detection threshold is initialized to 
0.5, then (SE/a) mtr = 0 and a detection occurs anytime the 
mean signal excess is equal to or greater than zero. In the 
probabilistic mode, the value ${SE/a} is compared to a random 
draw. If a -*■ 0, then for any non-negative value of SE, 

${SE/a} -*■ 1, and detections occur with probability 1. For 
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any negative value of SE, ${SS/a} -*■ 0, and detections occur 
with probability 0. Consequently, in the limiting case where 
a = 0, detections occur whenever the mean signal excess is 
equal to or greater than zero, and in this limit, the 
probabilistic mode is equivalent to the deterministic mode 
with the probability of detection threshold set tc 0.5. This 
suggests that the deterministic mode can be replaced by this 
limiting form of the probabilistic mode unless thexe is some 
utility in having (Pd) TH set to a value different from 0.5. 

The probabilistic mode corresponds to a complete inde- 
pendence model. In this mode, detection on a look is inde- 
pendent of detection on any other look. The cumulative 
probability of detecting a target in N looks is g;.ven by: 

N 

P{detecting in N looks} = 1-11 (l-Pd(i)} 1 (15) 

i=l 

-h 

where PD(i) is the detection probability on the i" look. 

For any non-zero probability of detection Pd, the cumulative 
probability of detection approaches one as the number of 
looks increase. Figure 7 shows the cumulative probability 
of detection for a complete independence model as a function 
of the number of looks, where the mean signal excess is 
constant and the single look probability of detection is 
0 . 1 . 

This figure also indicates the sensitivity of this model 
to the time between looks. It shows that after 22 looks, 
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Figure 7. Cumulative Probability of Detection Complete Independence Model 



the cumulative probability of detection exceeds 90 percent. 

It can also be shown that the average number of looks to a 
detection is 10. For this model to be valid, the time be- 
tween looks must be long enough to have the independence 
assumption hold. In the case of passive acoustic detections, 
the time between looks must also be long enough to satisfy 
the integration time requirements for detection threshold. 

In the NWCGS, the time between looks is a director adjusta- 
ble parameter initialized to 20 minutes. It is important 
that the game director understand how changing this value 
affects the cumulative probability of detection. The dis- 
cussion which follows describes some alternative models. 

1 . (X-q) Jump Model 

A more complicated model for the random component of 
signal excess than the one used in the NWCGS was introduced 
by J.D. Kettelle in 1960. it is referred to as the (X-a) 
jump model. In this model, the initial value of X(t) is 
determined in the same manner as that used in the NWCGS 
model; that is by a draw from a normal distribution having a 
mean of zero and a standard deviation of a . On subsequent 
looks, the value of X(t) remains the same until a jump occurs. 
The time of the jump is determined by a draw from an exponen- 
tial distribution having a mean of X. Each time a jump occurs, 
a new value for X(t) and a time to the next jump are deter- 
mined. The choice of lamda (known as the inverse relaxation 
constant) determines the model's level of dependence. As 
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A -*■ ", the mean time between jumps approaches zero, and the 
model approaches the independence model. As A -*• 0, the time 
between jumps approaches », and this limiting case corresponds 
to a complete dependence model. Complete dependence refers 
to the fact that if a detection does not occur on the detection 
opportunity, or ''look" with the largest mean signal excess, 
then detection will not occur on any other look. Therefore, 
the cumulative probability of detection for a series of looks 
is equal to the probability of detection for the look with 
the largest mean signal excess. Figure 8 shows a number of 
possible cumulative probability of detection curves available 
for a given mean probability of detection (0.1) and different 
values for A. 

2 . Gauss-Mar kov Models 

An alternative to the (A-a) jump model is the Gauss- 

Markov model. This model is Gaussian in nature because the 

joint distribution of a set of looks [X(l) ,X(2) , . . . ,X(n) ] is 

a multivariate normal distribution. For any X(i), its 

marginal distribution is normal with a mean of zero, and a 
2 

variance of c . It is Markovian, implying that given values 

for X(0) , X ( 1 ) , X(t-l), the value of X(t) depends only 

upon X ( t — 1 ) . This conditional dependence of X(t) on X(t-l) 

is normally distributed with a mean of p*X(t-l) , and a variance 
2 2 

of a *(l-p ). The parameter p, known as the autocorrelation 
coefficient, is defined by the equation: 

p = e~ At (15) 
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Figure 8. -Cumulative Probability of Detection Curves — Lamda-Sigma Jump Model 



The parameter A of equation (15) corresponds to the inverse 
relaxation coefficient of the (A-a) model. As A + 0 , the 
value cf p 1 and, as with the (A-a) model, the Gauss -Markov 
model approaches a complete dependence model. Similarly, 
as A -+■ 00 , the value of p -*■ 0 and the Gauss-Markov model 
approaches the independence model . The degree of dependence 
is determined by choice of A. 

Although the Gauss-Markov model is difficult to deal 
with analytically, a computer model that uses the conditional 
distribution of X(t) given X(t-l) is not difficult to develop. 
Cumulative probability of detection curves derived through 
simulation using various values of lamda are presented in 
Figure 9 . 

3 . K Out of N Models 

The k out of N detection model that is described here 
is based on independence model, however, each look results 
in either a success or failure rather than a detection or a 
ncn-detection. For a detection to occur, k successes must 
have occurred within the last N looks. Such models provide 
a means of describing operator characteristics. Figure 10 
shows how the k out of N model compares with the dependence 
and independence models of the NWCGS . The k out of N model 
offers the advantage of allowing one to set the time between 
looks to any specified time interval and then, by choice of 
k and N, of fitting the model to an appropriate cumulative 
probability of detection curve. Most detection processes 



34 




35 



Figure 9. Cumulative Probability of Detection Curves--Gauss-Markov Models 
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Figure 10. Cumulative Probability of Detection Curves — k Out of N Models 



require an initial integration time before a detection can 
be called, and this model can easily accommodate such a 
phenomenon. A k out of N model could also have been based 
on either the (A-a) jump model or the Gauss-Markov model. 

An example of a more general k out of N model can be 
found in the APSURF Surface Ship Engagement Model [Ref. 6]. 

In that model, active sonars were modelled as follows: 

1. A 6 out of 28 model determined detections if the 
operator was unalerted. 

2. A 5 out of 20 model determined detections if the 
operator was alerted. 

3. A 2 out of 10 model was used after a detection to 
determine if contact was still being maintained. 

This model is easily implemented and would require very 
little in the way of software modification to the existing 
NWCGS model . 

4 . Summary 

The alternative models described above are a sample 
of the models that have been used to describe the detection 
process. Although none of the models are without shortcomings, 
each of them has desirable properties not found in the current 
NWCGS models. In particular, they provide a range of serial 
correlation properties. This is appealing in modeling a 
detection process where the probability of detection time t 
may be influenced by information obtained at earlier times. 
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APPENDIX A 



The most effective mechanism for locating targets under 
the ocean surface is the detection and evaluation of acous- 
tic information specific to a target of interest. With 
this mechanism, detection may occur as a result of passively 
listening to the sounds generated by the target's propulsion 
equipment and auxiliaries, or may involve actively radiating 
acoustic energy in the suspected direction of the target and 
listening for the echo which would reflect off the target 
and return to the searcher. Passive sonar systems offer the 
advantage of being able to detect noisy targets at extremely 
long ranges (often in excess of one hundred miles) without 
alerting the target to the sdarch. Although active sonar sys- 
tems, or echo ranging systems, alert the target to what the 
searcher is about, they can be most effective against a 
quiet, slow-moving target close to the searcher. The system 
of choice is determined by the tactical situation and the 
environmental conditions under consideration. 

The passive acoustic modules used in the NWCGS are based 
on the generalized expected value sonar equations first for- 
mulated during the Second World War. Signal Excess, or the 
difference between the existing signal-to-noise ratio and 
that required to achieve a specified level of performance 
(Ref. 7] is the fundamental parameter that determines whether 
or not a detection has occurred in these modules. The required 



38 



signal-to-noise ratio (in decibels) is called the detection 
threshold. It is defined to be the signal-to-noise ratio 
that will give a specified probability of detection for a 
specified false alarm probability in a forced choice decision 
situation. If the specified probability of detection is equal 
to 0.5, detection threshold is usually referred to as recog- 
nition differential. 

After determining the appropriate detection threshold, 
the detection model determines the current signal excess 
level using the formula 



SE = SL - TL - NL + DI - DT 



(Al) 



for modeling passive detections, and 



SE = SL - 2TL + TS - NL + DI - DT (A2) 



for modeling active detections, where the parameters of the 
equations are defined as: 



SE - Signal Excess 

SL - The source level of the target as measured 
one yard from the acoustic center of the 
source 

DI - The directivity index of the receiver 

TL - The one-way transmission loss from the target 
to the source 

NL - The non-target noise existing in the area of 
the receiver 
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TS - The target strength as measured one yard from 
the acoustic center of the rarget 

DT - The detection threshold, or the minimum signal- 
to-noise ratio (in decibels) that is required 
to call a detection 

A more detailed discussion of these parameters follows. 

A. SOURCE LEVEL 

The value of the acoustic energy radiated into the medium 
by the searcher is termed the projector source level. For 
the sake of standardization, the source level is taken to be 
the acoustic power as extrapolated back to one yard from 
the acoustic center of the sound source. The interference 
patterns and other physical characteristics of acoustic 

energy close to its source produce measurements which are 

♦ 

of little value, so in practice the acoustic measurements 
are made from large distances, and the results can be 
reliably extrapolated to give reasonable source levels for 
all meaningful ranges. In passive acoustic listening, the 
signals of interest are those which are generated by the 
target alone. 

B. DIRECTIVITY INDEX 

Total acoustic power arriving at the receiving unit (be 
it a single hydrophone or array of hydrophones) consists of 
both signal power and noise power. Signal power is generally 
assumed to arrive at the receiver from a signal direction. 
Conversely, noise power is generally assumed to arrive at 
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the receiver uniformly from ail directions. A receiving unit 
that only accepts acoustic power from specific directions 
usually has a higher signal-to-ncise ratio than a receiving 
unit that accepts acoustic power omni-directionally. The 
signal-to-noise ratio gain due to the directivity of the re- 
ceiving unit is termed the directivity index DI of the 
receiver . 

C. TRANSMISSION LOSS 

As with any other form of radiated energy, acoustic 
signals are subjected to a number of physical effects which 
tend to distort and dissipate the signal. In the acoustic 
sonar models, these effects are collectively described by 
the parameter transmission loss (TL) . The primary constitu- 
ents of transmission loss are spreading, and absorption. 

An omni-directional sound source in a lossless infinite 
medium will radiate the same acoustic power over an ever- 
increasing area. Acoustic intensity, the quantity being 
measured by a detection system, is equal to power "P" divided 
by the total area "A" over which the intensity is distributed 

P = I * A (A5) 

If acoustic energy is subject only to spherical spreading in 
a lossless medium, the transmission loss can be described by: 

TL s = 20 LOG (r) (A6) 
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where r is the radius from the center of the sound source. 



In this case, intensity is inversely proportional to the 
square of the distance between detector and the center of 
the sound source . 

When physical or environmental boundaries serve to inhibit 
spherical spreading by restricting sound between two parallel 
planes (the ocean floor and surface as an example) , the power 
is radiated more like a cylinder of constant height "h" and 
increasing radius M r". The area of this cylindrical acoustic 
wave would be: 



Intensity is inversely preport. ional to the first power of 
range making the spreading loss component of transmission 
loss equal to 



Initially, the intensity of sound sources spreads spherically. 
At some distance r^_ away from the sound source, the acoustic 
intensity fills the medium between its bounds and spreading 
losses become more cylindrical than spherical. If the 
searcher is further away than this transition range, the 
transmission loss equation can be described by 



A 




(A7) 



TL 



s 



10 * LOG (r) 



(A8) 



TL 



10 * (LOG (R t ) + 10 * LOG (r) 



(A9) 
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The most common instance of this type phenomena both the 
target and the searcher are in the surface layer where the 
sound velocity profile is such that a lower acoustic "boun- 
dary" is formed. A similar, though less common example is 
the SOFAR, or deep sound channel formed when the target and 
searcher detection system are both located in a channel 
where sound is trapped from above and below by acoustic 
"boundaries" [Ref. 8], In either event, the threshold range 
R t is related to the height of the channel and depth of 
target below the top of the channel D t by 

R t = 1.094 * sqrt (D^/^-D^) (A10) 

where all quantities are measured in yards. 

1. Con v ergence Zone 

The nature of the oceans are such that given suffi- 
cient depth, all acoustic signals will be refracted back 
toward the surface. In this refraction process, the sound 
waves tend to converge, thereby reducing the losses due to 
spreading. These convergence rings, or zones as they are 
called, have transmission losses at least 10 dB less than 
they would be if the transmission were only spherical spread- 
ing. To account for this, transmission losses are calculated, 
and a convergence zone gain parameter CZ is subtracted from 
the expected transmission loss. 
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2 . 



Bottom Bounce 



The mud-sand composition of much of the ocean floor 
tends to absorb rather than reflect acoustic signals. If 
the bottom consists more of rocks and sand than mud, the 
reflective properties of those materials will allow trans- 
mission of sound at greater than anticipated ranges. To 
account for this effect, a bottom bounce adjustment factor 
BB can be subtracted from the transmission loss calculation. 

3 . Absorption Losses 

After spreading losses, the second greatest cause of 
intensity dissipation is absorption. This is caused by the 
conversion of acoustic energy into another form of energy. 

This type loss is usually attributed to shear viscosity, a 
friction-like conversion of acoustic energy into heat, or to 
energy involved in the dissociation and reassociation of 
ions in the ocean. The shear viscosity is associated with 
all frequencies of sound, while the dissociation reassocia- 
tion phenomena predominate at frequencies below 100 kHz. 

The intensity if absorbed exponentially with distance travelled 
through the absorbing medium and transmission losses asso- 
ciated with absorption can be described by: 

TL = ar (All) 

3 . 

where a is the absorption coefficient of the medium in units 
of decibels per kiloyard. A number of studies have been 
done in an effort to define an appropriate value for alpha. 
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It has been shown that the [Ref. 9] absorption coefficient 
is primarily a function of frequency and can be approximated 
using the equation: 

-2 SF T f2 -2 f 2 

a= 1.86 x 10 ^ + 2.86 x 10 z §- dB/Kyd (A12) 

+ f 2 T 

where f is the frequency of the sound, S is the salinity 
(in parts per thousand) of the sea, and F T is a relaxation 
frequency related to temperature T by the equation: 

F t = 21 . 9 . 10 (6-[1520/(T + 273)]) (a13 ) 

Additional causes of transmission losses are absorption of 
sourd into the ocean floor, scattering of sound energy at 

t 

the surface, leakage of sound energy out of a sound channel, 
changes in the absorption coefficient due to depth and tem- 
perature changes, etc. A complete accounting for all the 
elements which comprise transmission losses would be compu- 
tationally impractical and unnecessary for the model being 
discussed. As a minimum however, losses due to spreading 
and absorption should be considered. 

Transmission loss models are primarily a function of 
frequency of the sound source, and distance from the source 
to the target. As such, the one-way transmission loss from 
target to detector will be the same for active and passive 
sonars (assuming both systems are using the same frequency) . 
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Active sonar signals must travel from source to target and 
back, so the transmission loss in that circumstance is 
assumed to be exactly twice that of a passive sonar signal. 

D. NOISE LEVEL 

Any sound not originating from the target of interest 
is considered to be noise in the sonar model. Noise levels 
are generally categorized as either ambient noise (AN) or 
searcher self-noise (SN) . 

Ambient noise is the inherent noise of the ocean itself 
It is a collection of acoustic signals emanating from numer 
ous sound sources randomly arriving at the location of 
the searcher's hydrophone. It is primarily a result of the 
following factors: 

1. Environmental effects — tidal activity, rain, wind, 
seismic activity, turbulence, sea state, etc. 

2. Biologies — whales, dolphins, snapping shrimp, 
croakers, etc. 

3. Man-made effects — shipping intensity, petroleum 
exploration, deep water mining, etc. 

The diversity of sound sources involved in the makeup of 
ambient noise means that the entire acoustic spectrum is 
affected. It also implies that ambient noise is likely to 
be one of the most variable factors in the sonar equation. 
Ambient noise is usually modeled as a normal random varia- 
ble whose mean is dependent upon the shipping intensity in 
the area of the searcher and target [Ref. 10] . 
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Another major contributor to the undesirable noise at 
the detector is the self noise of the searcher. Auxiliary 
equipment is always being operated, and propulsion equip- 
ment is nearly always in use. As the speed of the searcher 
increases, so does the self noise problem associated with 
the propulsion equipment. At higher speed, the turbulence 
created by water flowing around the hydrophones becomes an 
additional source of self noise. Hydrophones are generally 
placed well forward of machinery spaces and have designed 
directivity away from any self noise sources, however total 
elimination of self noise is impossible. Self noise is 
usually modeled as a function of the speed of the searcher. 

E. TARGET STRENGTH 

I 

When a target is ensonified by an active detection sys- 
tem, a portion of the incoming power is reflected back to 
the searcher. If the target were a perfectly reflecting 
spherical body, the sound energy would be reflected uni- 
formly about the target, spreading spherically. In practice, 
targets are neither spherical, nor perfectly reflective. 
Target strength (TS) is the parameter that accounts for both 
absorption of energy by the target, and its reflectivity in 
the direction of the searcher. Target strength is strongly 
influenced by the relative angle of the target and searcher 
(aspect angle) . The variance of target strength with aspect 
differs widely among targets, but as a first order approxi- 
mation, Urick [Ref. 11] proposed an average submarine 
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target strength to aspect relationships which can be des- 
cribed by the equation: 

TS = TS *(16.17-2.98 cos (2B)-3. 083 cos (66) (A14) 

o 

where TS q is the maximum value of target strength, and 8 
is the aspect angle between target and searcher. Figure 11 
shows how this equation relates target strength to aspect 
angle . 

F . DETECTION THRESHOLD 

Target information available to the searcher is processed 

in some manner to arrive at a detection decision. The signal 

excess models used in the NWCGS assume that detections occur 

1 

whenever the signal-to-noise ratio (in decibels) exceeds some 
predetermined detection threshold DT. The choice of DT 
determines the probability of making a decision error. In 
a situation where a system is forced to make a decision after 
a single look, a detection will be called if the input exceed 
a threshold. When no target is present, the input is due 
to noise alone. In this event, if the input level exceeds 
the threshold, a detection will be erroneously called, and a 
decision error known as false alarm occurs. When a target 
is present, the input level is both signal and noise. In 
this event, if the combined signal and noise input level 
exceeds the threshold, the system will correctly decide that 
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Figure 11. Aspect Effect on Target Strength 



a target exists. Specifying the acceptable probability of 
false alarm determines a threshold. Having established a 
threshold, a specified probability of detection determines 
a required signal-to-noise ratio. This signal-tc-noise 
ratio (in decibels) is the detection threshold DT. 



50 



LIST OF REFERENCES 



1. Naval Postgraduate School Report NPS-51-79-006 , An 
Introduction to the Sonar Equations with Applications , 
by Alan B. Coppens , H.A. Dahl, and J.V. Sanders, pp . 35- 
50, February 1979. 

2. Naval Operations Analysis , 2nd ea. , p. 178, Naval Insti- 
tute Press, 1977. 

3. Sheehy, M.J. and Halley, H., "Measurement of the Attenua- 
tion of Low-Frequency Underwater Sound," Journal of 

the Acoustical Society of America , Vol. 29, pp. 464-469, 

4 April 1957. 

4. Urick, Robert J., Principles of Underwater Sound , 

pp. 96-104, McGraw-Hill Book Company, New York, 1975. 

5. Naval Postgraduate School Report NPS-61-79-006 , An 
Introduction to the Sonar Equations with Applications , by 
Alan B. Coppens, H.A. Dahl, and J.V, Sanders, p. 44, 
February 1979. 

6. 'Department of the Navy, Systems Analysis Office Report 

69-16, ASW Programs Surface Ship Engagement Model Abstract , 
by John D. Kettelle Corporation, p ." 12, January 1970. 

7. Naval Postgraduate School Report NP55Fo7i>041 , Some Notes 
on Search, Detection, and Localization Modeling , by 
Forrest, R. Neagle, pp. 4-13, April 1975. 

8. Cox, Albert W., Sonar and Underwater Sound , pp . 35-40, 
Lexington Books, D.C. Heath and Company, 1974. 

9. Urick, Robert J., Principles of Underwater Sound , pp. 96- 
103, McGraw-Hill Book Company, New York, 1975. 

10. NAVAIR Report 50-1G-24, ASW Oceanographic Environmental 
Services , V. 1, pp. 20.1-20.8, Naval Weather Service 
Command, Washington, D.C. 

11. Urick, Robert J., Principles of Underwater Sound , 

pp. 281-284, McGraw-Hill Book Company, New York, 1975. 



51 



INITIAL DISTRIBUTION LIST 



No. 



1. Defense Technical Information Center 
Cameron Station 

Alexandria, Virginia 22 314 

2. Library, Code 0142 
Naval Postgraduate School 
Monterey, California 93340 

3. Department Chairman, Code 55 
Department of Operations Research 
Naval Postgraduate School 
Monterey, California 93940 

4. Director of Wargaming 
Naval War College 
Newport, Rhode Island 02840 

5. Professor Neagle Forrest, Code 55Fo 
Department of Operations Research 
Naval Postgraduate School 
Monterey, California 33940 

6. Professor Michael Soverign, Code 55Zo 
Department of Operations Research. 

Naval Postgraduate School 
Monterey, California 9 3940 

7. Associate Professor Ray Shudde, Code 55Su 
Department of Operations Research 

Naval Postgraduate School 
Monterey, California 93940 

8. Associate Professor James Eagle, Code 55Er 
Department of Operations Research 

Naval Postgraduate School 
Monterey, California 93940 

7. LT David Martin Cashbaugh, USN 
USS John F. Kennedy (CV-67) 

Fleet Post Office 

New York, Mew York 09538 



Copies 

2 

2 

1 

1 

1 

2 

1 

1 

1 



52 




Thesis 
C27 475 
c.l 



20065] 

Cashba ugh 



if f L ' G 54 

?$ ocr P( 

~ U 1 V* • ' r* 



?8 ocr a 8 

:? r f 3 s; 



395 00 



f 9500 
/ 2 I i 2 



00 

02 2 



Thesis 
C27475 
c. 1 



200651 

Cashbaugh 

Examination of the 
sonar detection models 
used by the Naval War 
College gaming system. 



